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CMB polarization signal may be decomposed into gradient-like (E) and curl-like (B) mode. We have investigated E/B decomposition 
in pixel space. We find E/B mixing due to incomplete sky is localized in pixel-space, and negligible in the regions far away from 
the masked area. By estimating the expected local leakage power, we have diagnosed ambiguous pixels. Our criteria for ambiguous 
pixels (i.e. r^) is associated with the tensor-to-scalar ratio of B mode power spectrum, which the leakage power is comparable to. 
By setting to a lower value, we may reduce leakage level, but reduce sky fraction at the same time. Therefore, we have solved 
dACi/dr^ = 0, and obtained the optimal which minimizes the estimation uncertainty, given a foreground mask and noise level. We 
have applied our method to a simulated map blocked by a foreground (diffuse + point source) mask. Our simulation shows leakage 
power is smaller than primordial (i.e. unlensed) B mode power spectrum of tensor-to-scalar ratio r ~ 1 x 10"^ at wide range of 
multipoles (50 < / < 2000), while allowing us to retain sky fraction ~ 0.48. 
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1. Introduction 

Over the past years, CMB polarization has been measured by 
severa l experiments and is being measured by the Planck sur- 



veyor (' Kovac et al.ll2002l: iLeitch et alj|2002l: 



Halverso ~al.l 120021: iPrvke et alJ l2002bt 



^Vke et al."'2002a'; 
Leitchetal. 2005; 
Ade & et al. 2008; Prvke & et al. 2009; Hinderks & et al. 2008; 
iBrown & et al. 2009; The Planck Collaboration,2006.) . CMB po- 
larization pattern may be consider ed as the sum of gradient- 
like E mode and curl-like B mode (■Zaldarriaga & Seljay[1997l 
jKamionkowski et al. 1997). In the standard model, B mode po- 
larization is not produced by scalar perturbation, but solely by 
tensor perturbation. Therefore, measurement of B mode po- 
larization makes it possible to probe the universe on the en- 
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ergy scale at inflationary period (Kamionkowski et al. 19971 

jZaldai-riaga & Seliak 1997; Dodelson 2003; Liddle & Ly^ 
I2OO O; Mukhanov 2005). In most inflationary models, tensor-to- 
scalar ratio r is much smaller than one, and the WMAP 7 year 
data imposes an upper b ound on r < 0.36 a t 95% confidence 
level ( Larson et al.ll2010l:lKomatsu et al.ll2010l) . 

Besides instrument noise, there are complications, which 
limits detectability of tensor perturbation. Imperfection in re- 
moving foreground and gravitational lensing imposes observa- 
tional limit o n tensor-scalar-ratio: r ~ IQ~* and r ~ 3 x 10"^ 
respectively dSehak & Hiratal 120041: ITucci et al.ll2005h . Due to 
the nature of the observation or heavy foreground contamina- 
tion, reliable of estimation on CMB polarization signal is not 
available over a whole sky. Incomplete sky coverage leads to 
E/B mixing, and very signifi cantly limit our ca pacity to measure 
tensor perturbation as well (Bu nn et ani2003h . Therefore, there 
have been vari ous efforts to understand and reduce E/B mixing 
( Kim 200711 iBunn et aUlIool iLewis et alJl200l iLewisllIoof 
LSmith i20()6Y 

It is best to implement E/B decomposition in map space, 
since diffuse foregrounds and point sources are well-localized 
in map space, and their spatial information are known relatively 



better than other properties. In this paper, we investigate E/B 
decomposition in pixel space. Our investigation shows that E/B 
mixing is highly localized in pixel space. Therefore, we may re- 
duce E/B mixing effectively by excluding the ambiguous pixels. 
We have applied our method to simulated maps partially blocked 
by a foreground (diffuse + point source) mask. After excluding 
ambiguous pixels, we find that leakage power in retained pixels 
(sky fraction ~ 0.48) is smaller than primordial (i.e. unlensed) B 
mode power spectrum of tensor-to-scalar ratio r ~ 1 x lO"-' at 
wide range of multipoles (50 $ / $ 2000). 

The outline of this paper is as follows. In Sec.|2l we discuss 
all-sky analysis of CMB polarization. In Sec. [3] we derive E/B 
decomposition in pixel space. In Sec.lH we discuss the applica- 
tion to cut sky, and the method to diagnose ambiguous pixels. 
In Sec. |5]and|6] we present our simulation result. In Section |7] 
we summarize our investigation. In Appendix|Al we discuss er- 
ror analysis of pseudo C/ estimation, and show interpixel noise 
correlation may be neglected. 



2. STOKES PARAMETERS 

The state of polarization is described by Stokes parameter 
(Kraus 1986; Rohlfs & Wilson 2003). Since Thompson scat- 
tering does not generate circular polarization. Stokes param- 
eter Q and U are sufficient to describe CMB polarization 
dPodelsonl l2003h . Stokes parameter Q and U transform under 
rotation of an angle t/* on the plane perpendicular to direction n 
(Zal dan-ia ga & Seljak 1997; Za ldarriaga 1998.) : 



(Q + lUYin) = e^^"^{Q + iU)(n). 



(1) 



Therefore, all-sky Stokes parameters may be decomposed into 
spin ±2 spherical harmonics (Zaldarriaga & Seliak,1997.) as fol- 
lows: 
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Q(n) + iU(n) = ^ fl±2,/m ±2J'/m(n), 



(2) 
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where the decomposition coefficients fl±2,/m are obtained by: 

«±2,/„, - j[Q(n)±iU(h)]^2Y;jA)dn. (3) 

Though the quantity shown in Eq. [2]has direct association with 
physical observables (i.e. Stokes parameters), rotational variance 
leads to computational complication. Therefore, two real scalar 
quanti ties, termed 'E' and 'B' mode, are often built out of Q(n ) + 
iU(n) dKamionkowski et al.|[T997tlZaldarriaga & Seliaklll997l) : 

£(n) = -l[SHQ(n) + iUm + S\Q(n)-iUm], 



Im 



(1 + 2)1 
\l-2y. 



(4) 



B(fl) = -[S\Qin) + iUin))-d\Q{n)-iU{m, 




(5) 



where d and d refer to lowering and raising operator respec- 
tively (IZaldarriaga & Seliak 1997h. The explicit express ion of d 
and S are given as follows ( Zaldarriaga & Seliaklll997l) : 



■ sm' 



6 j(e,^) = -sin^'e 



d d 
— + I CSC 9— 
86 d(p 

' a d 
I CSC ti- 
de d(p 



sin-' j{e. 



sin' 6 J(0, 



where (p) is an arbitrary spin s function. The de- 

comp osition coefficients of E and B mode are related to 
a4.7 /m dZaldarriaga & Seliak|[l997h as follows: 



fl£,/m - -(fl2,/m + fl-2,;m)/2. 



(6) 
(7) 



For a Gaussian seed fluctuation model, decomposition coef- 
ficients of E and B mode satisfy the following statistical proper- 
ties: 



(<^*B.lm"B,fm'} 



; Oii'On 



where (. . .) denotes an ensemble average. In Fig.[T] we show un- 
lensed Cf ^ and Cf * of the WMAP concordance ACDM model 
for various tensor-to-scalar ratio r. 



3. E/B decomposition in pixel space 

In this section, we are going to derive a pixel-space analogue of 
E/B decomposition. Using Eq. [3]|6]and|7] we may easily show 
Eq.|4]and|5]are equivalently given by: 



E(n) 




(1 + 2)1 
(I -2)1 



CEJm Ylm(n) 



(l + 2)\ 



(1-2) 



I "EJm Im 



YL(n), 



= -^(/ F^(h',h)[Q(n')-iU(n')] dQ' 
+ J F^(n', n) [Q{n') + iU(n')] dQ' 



(10) 




1000 



1000 



Fig, 1. the power spectrum of E (top) and B (bottom): no lensing, 
B mode power spectrum is plotted for various tensor-to-scalar 
ratio r. 




BJm Im 



= ^JJ F,(n',n)[Q{h')-iU(h')]dQ' 
- j F4n', n) [e(n') + iU(h')] dQ'j , 



where 
F^(n',h) 



Im 



(1 + 2)1 
(1-2)1 



±2Y,„,(h') F* (n) 



(11) 



(12) 



Therefore, we may identify F±{h',n) as pixel-space filters for 
E/B decomposition. Using the property of spin-i spherical har- 
monics sY^^^(h) = -sYi-m(fi), we may show the pair of the filter 
functions have complex conjugate relation: 



F+(n',n) = F*_(h',n). 



(13) 



In Fig. |2] we show the real and imaginary part of F+{n', n) for 
a fixed n. In Fig. |3] we show one dimensional plot of \F±{n', n)| 
for a fixed n. As shown in Fig.|2]and|3] our pixel-space filter pos- 
sesses sharp peaks around n. Note that Eq.[T2] would approach 
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i(02,lm - «-2,;m)/2, 




Fig. 2. Filter function: Re[F+(n',n)] (left), Im[F+(n',n)] (right) 
for a fixed n, and n' spanning 2° x 2° area, F+{n',n) - 
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Fig. 3. Filter function: modulus |F±(n',n)| for (6 = 7r/2,</> = 0), 
highly peaked at cos 6' = 



6{n' - n), if ±2i'/m(n') were Yi,„(n'). Using Eq.\U\ we 

may easily show Eq.[TO]and[TT]are equivalently given by: 

£(n) = - J dn'Re[F4h',A)iQ(h')-iU(h'))], 
= - J dQ'Re[F_(n',fi)(e(fi') + /t/(n'))]- 

B(A) = - J dQ'Im[F+(n',n)(e(fi')-/f/(ft'))], 
= J" dQ'Im[F_(n',n)(e(n') + /f/(n'))]- 

4. Incomplete sky coverage 

Due to heavy foreground contamination, CMB polarization sig- 
nal is not estimated reliably over a whole sky. For instance, the 
WMAP team have subtracted diffuse foregrounds by template- 
fitting, and masked the regions that cannot be cleaned reliably. In 
Fig. 121 we show a foreground mask, which combines the WMAP 
team's polarization mask with the point source mask and shall 
be used for our simulation. The E/B decomposition coefficients 
from a masked sky are given by: 



fl£,/m - -(aijm + S-2,;m)/2, 



(14) 



where 



= J Win' 



)[Q(h)±iU(h)] ^2YL{n)d^l. 



(15) 



(16) 



and W(n') is a foreground mask. Therefore, E/B maps recon- 
structed from incomplete sky are given by: 



E(n) = 



(1 + 2)1 



V(/-2)! 



(17) 



-i J J W(n') F,(n', h) [Q(n') - iU(n')] dQ.' 
J Win') F.in', n) [Q(h') + iU(n')] dfi' 



B(n) = 



A/(7^"*'""^""^"^' 



(18) 



= ^ J J W(n') F^(h', n) [Q(h') - iU(n')] dQ' 
- J W(h')F4n',h)[Q(h') + iU(h')]dQ' 



Since filter functions F±(n', n) are sharply peaked around n, cer- 
tain pixels far away from masked regions may contain negligible 
E/B mixing, and vice versa. Equivalently, E/B mixing is local- 
ized in pixels close to the masked regions. For higher pixel res- 
olution, Eq.[T2]contains summation up to higher /, which makes 
the peak of the filter function sharper. Therefore, E/B mixing de- 
creases with increase in pixel resolution. The spherical harmonic 
method (the first line of Eq.fTTjandfTsTl are much faster than the 
pixel-space method (the second line), while they are mathemat- 
ically equivalent. Therefore, we are going to rely on spherical 
harmonic transformation method for our simulation in the next 
section. However, it should be kept in mind that the pixel-space 
approach have provided useful insights on E/B decomposition of 
incomplete sky. 

Using Eq. [18] we may show the expected power of B(h) is 
given by: 



(bHa)) 



dQ" W(n')W(n")x 



(19) 



(FJh',n)F4n",n){(Q(n') - iU(n'))(Q(n") + iU(n"))} 
+F4n',h)F^(h",h){(Q(h') + iU(h'))(Q(h") - iU(h"))} 
-F^(h',h)F4n",n){(Q(n') - iU(h'))(Q(h") - iU(n"))} 
-F_(n',n)F_(n",n)((e(n') + iU(h'))(Q(h") + iU(n")))). 

where (. . .) denotes an ensemble average, and 

<(e(n') + iU(n'))(Q(n") + iU(n"))) 



2/+ 1 

4-K 



(Cf^ + Cf") 2Yi,-2(Ii,0)e 



-^BB 



2i(a+y) 



{(Q(n') + iU(n'))(Q(n") ± iU(h"))} 

= Z 



(20) 



(21) 



2/+ 1 



4n 



(Cf^ + Cf) 2l',,-2(j8,0)e^2/(<.-r)^ 



and /3 is the separation angle between n' and n", a and y are 
the rotation angles respectively, which align eg at n' and «"with 
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Fig. 4. effective sky fraction /sky for various cut level r, 



the great c ircle p assing through h' and h" (refer to Fig. 1 of 
iNg & Liu ('1999) for a geometrical diagram). Taking into ac- 
count Eq. [191 [20] and [21] we may easily show that the expected 
leakage power at n is given by 



n)) ^IJdn' 



dO" W(n')W(h")x 



iF^(n',n)F4n",mQE(n') - iUE(fi')){QEin") + iUE(A"))} 
+F4h',n)F4h",mQE(n') + iUE(n')){QE(n") - iUE(A"))) 
-F^in',n)F^(n",mQE(n') - iUE(n'))(QE(n") - iUE(n"))) 
-F4n',n)F4n",mQEin') + iUEin'))(QEin") + /{/^(n")))). 



where 

{(Qe(A') + iUE(n'))(QE(n") + iUE(n"))) 

= Z 



(23) 



2/+ 1 
47r 



Cf^ 2Yi,-2(P,0)e 



,2i(ff+r) 



{{QEin') + iUEih '))(QE(n") ± /f/£(n"))) 



(24) 



Therefore, we may diagnose ambiguous pixels (i.e heavy E/B 
mixing) by comparing (B|(n)) with (B|(n)), where (Bg(n)) is 
the local power contributed by B mode and given by replacing 
'E' with 'B' in Eq.[22l[23]and[24| However, estimating Eq.l22l 
[23] and [24] is prohibitively complicated. Therefore, we are go- 
ing to resort to Monte-Carlo simulations in order to estimate 
(B|(n))/(B|(n)>. Depending <B|(n))/(B|(fl)), we may classify 
the pixel at n as 'pure' and 'ambiguous'. To be specific, we may 
retain pixels satisfying: 



(Blih)} ^ r 



(25) 



where r is the assumed tensor-to-scalar ratio of Monte-Carlo 
simulation, from which (B|(n))/(B^(n)) is estimated. Therefore, 
the level of leakage in retained pixels is comparable to the pri- 
mordial B mode power spectrum of tensor-to-scalar ratio r^. . In 
Fig. [H we show the sky fraction for various r^, given a fore- 
ground mask shown in Fig.[2l Since sky fraction decreases with 
lower Tc, we may not simply set r^. to a lowest value. Therefore, 
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Fig. 5. Numerical solution of Eq.[27|for various r and the noise 
level of Planck HFl instrument: two plots represent the Left 
Hand Side (LHS) and Right Hand Side (RHS) of Eq.[27l 

we need to derive an optimal rc, which minimizes the estimation 
error. The estimation error of B mode power spectrum is given 
by: 



(22) 



ACf" = 



-(c, 



(2/+l)/sky 
(2Z+l)/,ky' ' 



BB ^ gEE 



+ Ni), 



+ -C,+ Ni), 



(26) 



where A^; is noise power spectrum. Note that the leakage does 
not bias the B mode power spectrum estimation, but increases 
the variance, when the power spectrum estimation is made by 
a pseu do-Q method and leakage is taken care of (Hivon et aTj 
l2C)02aLlGrainet al.il2009h . By requiring 5 ACf^/Sr^ = 0, we get 



5 In/, 



sky 



r + r.+Ni/Cf^ir = 1)' 



(27) 



In Fig. [5] we plot the left and right hand side of Eq. [27]for the 
noise level of Planck HFl instrument, and the multipole / = 86, 
which is the peak multipole of primordial B mode power spec- 
trum. From Fig.[5] we find curves intersect at r^. a; 4 x 10"^ with 
weak dependence on r. It should be noted that the weak depen- 
dence is due to the low signal-to-noise ratio of the considered 
experiment (i.e. Ni/Cf^(r = 1) » 0), and the dependence on r 
is not weak in general. We are going to use ~ 4 x 10"^ for the 
simulation in the next section. 



5. Application to simulated data 

Using the WMAP concordance ACDM model, we have sim- 
ulated Stokes parameter Q and U over a whole-sky with a 
HEALPix pixel resolution (Nside=1024) and 10' FWHM beam. 
We have made the inputmap to contain no B mode polariza- 
tion. Therefore, any non-zero values in output B map are at- 
tributed to leakage. We show our simulated polarization map 
in Fig. [6] where the orientation and length of headless arrows 
indicates polarization angle and amplitude respectively. Note 
that the polarization map shows only gradient-like patterns, be- 
cause they contain only E mode polarization. It is well-known 
that E/B mixing i ncreases with the length of cut sky boundary 
dBunnet al.ll2003l) . We have combined the WMAP team's po- 
larization mask with the point source mask, and prograded it 
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Fig. 6. Input polarization map: E mode polarization only 




Fig. 7. foreground mask: smoothed with 1.5° FWHM Gaussian 
kernel. 




Fig. 8. output B map from incomplete sky 



to Nside=1024. In order to reduce sharp boundaries, we have 
smoothed the mask with 1.5° FWHM Gaussian kernel. We 
have referred to the WMAP team's bou ndary smoothing proces s 
of Internal Linear Combination map (iHinshaw & et alJl2007h . 
Nevertheless, it should be noted that smoothed boundary is not 
essential to our method, and further improvement m ay be pos- 
sible by using more sophisticated smoothing kernel dPas et alJ 
I2OO9I) . In Fig. 13 we show our smoothed mask, whose sky frac- 
tion amounts to 0.71. 

In Fig. [8] we show a B map, which we have produced from 
the masked polarization map. Using (B|(n))/(B|(n)) estimated 
from lO"' simulation, we have diagnosed ambiguous pixels, and 
retained pixels of <B|(n)>/(B|(n)) < 4 x IQ-^/r. In Fig.|9l we 
show the B map, where ambiguous pixels are excluded. We find 
retained pixels of Fig. |9] amount to sky fraction /s^y = 0.48. 
From the retained pixels, we h ave estimated the leakage power 
spectrum by pseudo C/ method dWandelt et al.l200"lllHivon et alJ 




Fig. 9. filerted B map: masked by a foreground mask, and pixels 
of <B|(n))/<B|(n)) > 4 x 10"^ are set to zero. 





leakage (no filtering) 




leakage (filtering) 




r=10"'' 




r=10"^ 









500 1000 1500 2000 
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Fig. 10. leakage power spectrum and primordial B mode power 
of various tensor-to-scalar ratio r. a blue curve correspond to the 
leakage power estimated without ambiguous pixel filtering (Fig. 
[8]l, a green curve to the leakage power estimated with ambiguous 
pixel filtering (Fig.|9]l. 

l2002bh . Power spectrum is usually estimated by a pseudo C/ 
method at high multipoles (/ > 30), while by maximum like- 
Uhood method or Gibbs sa mpling at low multipoles (I < 30) 
dBond et al. 1998; Eriksen eTaP 120041: iHinshaw & et al J 12007 : 
fefstathiou 2006). However, we find pseudo C/ method at low 
multipoles is good enough for our need, since we do not in- 
tend accurate estimation of likelihood function. Besides that, we 
are mainly interested in leakage power at (/ > 30), because pri- 
mordial B mode power spectrum has a peak around multipoles 
/ ~ 90 (see Fig.[TJ. In Fig. [10] we show the leakage power spec- 
trum and B mode power spectrum of various tensor-to-scalar ra- 
tio r. Fig. [To] shows leakage power (green) is smaller than B 
mode power spectrum of r = 10"^ at wide range of multipoles 
(50 ^ I ^ 2000), when ambiguous pixels are excluded (i.e. the B 
map in Fig.|9]l. In Fig. [10] we also show the leakage power (blue) 
estimated without filtering (i.e. the B map in Fig.|8]l. Obviously, 
we have reduced leakage significantly by excluding ambiguous 
pixels. 

6. scale-dependence of leakage 

In order to reduce leakage, we have removed ambiguous pix- 
els in a scale-independent way. On the other hand, it is known 
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Appendix A: error analysis 

Power spectrum is usually estimated b y pseudo-Q method 



at high multipol es ( Wandelt et alj 200ll : Hivon et al 



that leakage has some dependence on scales as well as real- 
space. Specifically, leakage of low I extends over large area, 
while leakage of high / is often confined to the small area nearby 
the boundary. Therefore, one may argue that our method does 
not reduce the leakage at low / as effectively as that of high I 
or leads to unnecessary loss of information. However, as shown 
in Fig. [8l leakage is relatively localized in pixel space, and we 

were able to reduced leakage very significantly at wide range _ 

of multipoles, while retaining sky fraction 0.48. Besides that, jEfstathioul 120061: iNolta & et alJ l2008t Eson et al 
our simulation shows the leakage at lowest / is reduced signifi- 
cantly as well. It is also possible to implement a further leakage 
reduction in a scale-dependent way, after ambiguous pixels are 

removed. Nevertheless, a hybrid method, which exploits scale ^bb _ ^^(m-I)^^,^^* 
and position dependence simultaneously, may be most optimal. ' 
A wavelet approach may be promising for such implementa- 
tion, since wavelet functions are, in general, well-localized in 
harmonic space as well as pixel space. We defer a rigorous in- 
vestigation to a separate publication. 



According to pseudo-C/ method, we may estimate power spec- 
trum as follows (Wand elt et al.,.2001;,Hivon et aL,2002b.) : 



where 



2Z' -H 1 , 

M„, = — ; 2j I^'" 



7. Discussion 

We have investigated E/B decomposition in pixel space, and 
shown that we may produce E/B decomposed maps by convolv- 
ing polarization maps with certain filter functions of a sharp 
peak. We find that E/B mixing due to incomplete sky is local- 
ized in pixel-space, and negligible in the regions far away from 
masked area. By estimating the expected local leakage power 
and comparing it with the expected pure mode power, we have 
diagnosed ambiguous pixels and excluded them. Our criteria for 
ambiguous pixels (i.e. r^) is associated with the tensor-to-scalar 
ratio of B mode power spectrum, which the leakage power is 
comparable to. The estimation error AC/ may increases with 
lower Tc, because sky fraction decreases. Therefore, we have 
solved dACi/drc - and obtained the optimal r^, which min- 
imizes the estimation error, given a foreground mask and noise 
level. We have applied our method to simulated maps blocked by 
a foreground mask. Simulation shows that leakage power is sub- 
dominant in comparison with unlensed B mode power spectrum 
of r ~ 1 X 10"^ at wide range of multipoles (50 ^ I ^ 2000), 
while pixels of sky fraction 0.48 are retained. We may apply our 
method equally to small sky patch observation, by treating unob- 
served sky as masked region. From simulation with sky patch of 
simple symmetric shape, we have confirmed our method reduce 
E/B mixing very effectively. A rigorous investigation is deferred 
to a separate publication. 

Noise is slightly correlated from pixel to pixel in E and B 
maps, even when interpixel correlation is absent in Q and U 
maps. However, this interpixel noise correlation induced by E/B 
decomposition is not confined to our method, but E/B decompo- 
sition in general. Besides that, we find interpixel noise correla- 
tion may be neglected without sacrificing the accuracy of error 
analysis (refer to AppendixlAlfor details). Therefore, it does not 
limit the applicability of our method. 

Current observations such as WMAP were unable to detect 
B mode polarization. Therefore, we did not attempt to apply our 
method to observation data. When Planck polarization data of 
high Signal-to-Noise-Ratio (SNR) are available in near future, 
we may apply our method to the data, and be able to detect B 
mode polarization. 
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1 I' I" 

2 0-2 



I I' I" 
-2 2 



and 



Wl» 



-f 



df2W(n)y*„„(n), 



(A.l) 



(A.2) 



(A.3) 



with W(n) being a foreground mask function. The pseudo- 
quantity Ci and a{ are given by: 



Ci = a! (, 

' 21+1 ^ 



air. 



Im J 



(A.4) 



(A.5) 



where j and k refers to a DA-year combination and A(n) refers 
to data (T, E or B). We may split a-j^^^ into signal and noise: 

al = a,,„+Nl^. (A.6) 
The noise part is given by: 

47T 



Wpix ,. 



(A.7) 



where a pixel index / runs over pixels outside a foreground 
mask, and A^,- refers to noise at ith pixel. By central Uniit the - 
orem (lArfken & Webed 120001: Ik. F. Rilev M. P. Hobsonll2006h . 
Nj^^^ follows a Gaussian distribution, and is uncorrelated among 

distinct ;th DA-year data (i.e. {NjjNiJ*} Sjk). 

If cross power spectra are used (i.e. j + k), noise does not 
bias estimation, and its stati stical pro perties need to be known 
only for errors analysis (Hi nshaw & et al. 2007). In order to un- 
derstand the effect of noise on error analysis, let us consider co- 
variance of C;: 



Cov(C,,,C,J = (Q,C,,)-Q,C,,, 

where C; is a theoretical power spectrum and 

(CM = (y(M-'),,,Qy(M-')y, c,>. 



(A.8) 



^(M-i),,,^(M-')y-<C,Q). 



(A.9) 
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Using Eq. lA.4l and lA.6l we find 

(C,C,)=A (A. 10) 



1 



(21+ 1)(2/' + 1) 

where A denotes terms irrelevant to noise. Therefore, we need to 
estimate noise covariance {Nl^iNj,^^,)*} in order to estimate Eq. 
lA.lOl First, let us consider diagonal elements of noise covari- 
ance: 



(A. 11) 



«2. 

pix ; 



pix ; i'>i 



The main contribution of Eq. lA.llI comes from the first term, 
because of cancellation through summation in the second term. 
Therefore, we find with good approximation: 



^ * 



pix I 



OfF-diagonal elements of noise covariance are given by: 



n ■ ^ 

pix I 



(A.12) 



pix / i'>i 



Comparing IA.12I with Eq. lA.l II we may see the magnitude of 
off-diagonal elements is much smaller than that of diagonal ele- 
ments, because cancellation through summation arise both in the 
first and the second term of Eq. IA.12I Therefore, we find noise 
covariance as follow: 



167r- 



(A. 13) 



pix I 



Using Eq. lAiSirOllAlOl and lAJJl we find covariance of Q: 
Cov(Q,,4) « C + 2(M-'),,,(M-i)w (A. 14) 



pix^ ^ m ii' 



where C denotes terms irrelevant to noise. As shown in Eq. lA.14| 
we may neglect interpixel noise correlation in computing covari- 
ance of Ci, 
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